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The introduction of metasurfaces has renewed the SnelPs law and opened up new degrees of freedom to 
tailor the optical wavefront at will. Here, we theoretically demonstrate that the generalized SnelPs law can be 
achieved for reflected acoustic waves based on ultrathin planar acoustic metasurfaces. The metasurfaces are 
constructed with eight units of a solid structure to provide discrete phase shifts covering the full In span with 
steps of ft/ 4 by coiling up the space. By careful selection of the phase profiles in the transverse direction of the 
metasurfaces, some fascinating wavefront engineering phenomena are demonstrated, such as anomalous 
reflections, conversion of propagating waves into surface waves, planar aberration-free lens and 
nondiffracting Bessel beam generated by planar acoustic axicon. Our results could open up a new avenue for 
acoustic wavefront engineering and manipulations. 
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The angles of reflected and refracted waves are certainly determined when light impinges on a planar interface 
between two media with different refractive index due to the conservation of momentum along the tan- 
gential directions of the boundary, as is well known from the famous Snell's Law. Recently, the Snell's law 
was revisited in the context of metasurfaces with phase discontinuities constructed by metallic antennas capable 
of providing discrete phase shifts covering 2n span 1 . The generalized Snell' law was introduced by applying 
Fermat's principle to interpret the anomalous reflection and refraction phenomena. Inspired by this pioneering 
work, fascinating optical wavefront engineering capabilities, such as light bending 2 , anomalous reflection 3 , con- 
verting propagating waves to surface waves 4 , flat lens and axicons 5 , etc., have been demonstrated both theoret- 
ically and experimentally through metasurfaces with phase discontinuities 6 . 

As another classical waves, acoustic waves also follow Snell' law. It is therefore apparent that the generalized 
Snell's law and metasurfaces, if could be successfully realized in acoustics, would have deep implications for 
acoustic devices, acoustic applications and the field of acoustics in general. A straightforward analogy between the 
optical and acoustic metasurfaces, however, is infeasible due to the inherent distinction between the electromag- 
netic and acoustic waves. In optics, the concept of metasurface is generally implemented by utilizing the plas- 
monic resonance mechanism through metallic antennas, which is able to tailor the wave phase from 0 to 2n but 
has no direct counterpart in acoustics so far. On the other hand, the use of plasmonic resonance inevitably leads to 
a direct result that the optical metasurface can hardly have a deep -sub wavelength size and operate only with cross- 
polarized fields 7 . It is therefore necessary to explore intrinsically new physical mechanism in an attempt to achieve 
the acoustic metasurfaces with ultrathin structures, which will be particularly significant for the practical applica- 
tions of the resulting devices. Considering the fact that the naturally available materials in acoustics can only 
provide limited phase shifts, the concept of generalized Snell's law cannot be analogously translated from 
electromagnetic waves to acoustic waves. The emergence of acoustic metamaterials 8 " 19 , which could provide 
unavailable parameters in nature, has significantly broadened the horizon for acoustic waves. Novel phenomena, 
such as extraordinary acoustic transmission 20 " 24 , sub -diffraction imaging 25 " 27 , cloaking 28 ' 29 , one-way transporta- 
tion 30 " 33 , super absorption 34 , etc., have been demonstrated both theoretically and experimentally. Hence, the 
advance in acoustic metamaterials may provide possibility and should be good candidates to realize the acoustic 
metasurfaces. 

In the following, we show that generalized Snell's law for acoustic reflected waves can be realized through 
ultrathin planar acoustic metasurfaces by coiling up space with two stiff corrugated beams. By appropriately 
selecting the geometrical parameters, eight units can support discrete phase shifts ranging from 0~27i with steps 
of n/4. Anomalous reflections, conversion of propagating waves to surface waves, flat acoustic aberration-free lens 
and nondiffracting Bessel beam generated by flat acoustic axicon are also demonstrated by imposing suitable 
phase profiles in the transverse direction of the metasurfaces. Our results should open up new degrees of freedom 
in acoustic wavefront engineering, such as the field of caustic engineering. 



SCIENTIFIC REPORTS | 3 : 2546 | DOI: 1 0.1 038/srep02546 



1 




0.45 0.50 0.55 0.60 

a [cm] 



Figure 1 | An acoustic metasurface for generalized Snell's law. (a) The schematic diagram of an acoustic metasurface made of two stiff corrugated 
beams with a channel coiling up the space. The coiling structure has a width p = 1 cm and length a(a = 0.8 cm in this example). The width of the channel 
is d = 0.067p. The width of the beams is w = 0.03p, and the corrugation length is / = a — 2 w — d. Sound hard boundary conditions are imposed to the left 
boundaries (red line) to mimic the fact that the metasurface is actually coated on a stiff plate. The light blue and dark red arrows refer to the propagation 
direction of incident and reflected waves, respectively. The label "1", "2" and "3" refer to the three outlets of each element, (b) The phase of the 
reflected waves, as a function of the length a of the metasurface, with incident wavelength X = 1 9.6p. The black dots refer to specific a values for eight units 
to fulfill the desired discrete phase shifts, (c) Schematic diagram of the eight units with the specific a values shown in (b) with black dots. The gap between 
each unit is p— {a\ + a 2 )/2, with a x and a 2 denoting the length of the two adjacent coiling structures, (d) The pressure strips of the reflected waves 
by the eight units. The high maps of pressure field are utilized to clearly show the different phase shifts by each unit. 



Results 

The acoustic metasurface, as shown in Fig. 1(a), is constructed by 
immersing identical thin rigid beams (width w and length /) in air in a 
way that forms fluid channels (width d). It has been proven that the 
coiling structure (width p and length a) could behave as an effective 
medium with extremely high refractive index 18,35 , such that the 
reflected waves can be delayed heavily after the incident waves 
impinge on the structure and propagate within the zigzag path. It 
is therefore reasonable to infer that the phase of the reflected waves 
could be controlled at will by appropriately selecting the geometrical 
parameters of the structure, and will be demonstrated later. In order 
to ensure a planar structure, the length a of the metasurfaces assoc- 
iated to the length / of the rigid corrugations are tailored to yield the 



desired phase shifts, while the values of w, d and p are fixed in the 
calculations. The phase of the reflected waves as a function of a with 
free space wavelength X= 19.6 cm is plotted in Fig. 1(b). It should be 
noted that, since acoustic waves, as scalar waves, can propagate 
within the channels freely, the reflected waves at the right-edge of 
the metasurface (containing three channels labeled with "1", "2" and 
"3" in Fig. 1(a)) can be considered as the accumulated field of three 
sources, resulting in the inhomogeneous pressure distribution at the 
right edge of the metasurface. Therefore, the phase information at the 
right edge of the metasurface is retrieved by integration along extra 
line X/8 (not shown in Fig. 1(a)) away from the edge. From Fig. 1(b), 
it is sufficient to frame eight units that could realize discrete phase 
shifts covering the full 2n span with steps of n/4. The exact values of a 
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Figure 2 | The generalized Snell's law of reflection and anomalous reflections, (a) Schematics for the derivation of the angle of reflection, (j) and 
(f> + d§ are the phases at the two cross points separated by dy along the y direction. 6 r represents the anomalous reflection angle induced by the discrete 
phase shifts, (b) Pressure field pattern for the gradient phase profile of d(f>/dy = n/\6p. The black arrows refer to the theoretical value of the 
reflected angle. 
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Figure 3 | Metasurfaces for converting propagating waves to surface 
waves. Pressure field pattern for the conversion of propagating waves to 
surface waves on the metasurface with phase profile (j)(y) = — 2ny/X. 
Surface waves can propagate near the metasurface with high -efficiency 
conversion. The black arrows indicate the wave vectors of the surface 



for achieving these discrete phase shifts are also illustrated with black 
dots in Fig. 1(b). Schematic diagrams of these eight units with dif- 
ferent values of a are shown in Fig. 1(c). Planar acoustic metasurfaces 
are obviously constructed because the lengths of the units in x dir- 
ection are entirely identical. To further verify the discrete phase 
shifts, the reflected waves by these eight units are shown in 
Fig. 1(d). The strips refer to the pressure filed patterns p(x,y) at the 
same time instant. The peak of the pressure field can shift up to a 
wavelength, which provide solid support that the phase shifts cover 
the whole 2n range. It therefore could be concluded that desirable 
discrete phase shifts can be realized by the eight units. The metasur- 
faces constructed by the eight units are ultrathin with the thickness 
approximately equal to X/19.6, revealing that the metasurface is good 
candidate to realize acoustic devices easy for integration. It is worth 
pointing out that fascinating features of the coiling metamaterials 
such as negative refractions and zero refractive index are first pre- 
sented theoretically and later verified experimentally 18 ' 36 ' 37 . Recently 
we have proved that by coiling up space it is possible to change the 
propagating phase of acoustic wave and thereby focus the waves, but 
the phase delay in that work is too limited to cover 2n range 35 . In what 
follows we will further demonstrate the potential of coil-up acoustic 
metamaterials to design a more general metasurface that has an 
ultrathin structure and is capable of controlling reflected acoustic 
waves arbitrarily. 

The realization of In range phase shifts for the reflected acoustic 
waves allows us to revisit the Snell's law by selecting appropriately 
phase profile. Figure 2(a) illustrates the schematic diagram for the 
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Figure 4 | Planar acoustic lens constructed by an acoustic metasurfaces. (a) Schematic diagram of the design of lens. A hyperboloidal phase profile 
along y direction is utilized to focus acoustic plane wave to a single point at a distance /from the metasurfaces. The green spherical line with radius /is the 
desired equiphase surface of the reflected wave so that the phase shift at P L should be proportional to the length of the red line PlS l . The theoretical 
continuous phase shifts (red dots) and the discrete phase shifts provided by the metasurfaces (blue squares) along the y direction are shown in (b). 
The half-height of the lens is chosen to be 80 cm. (c) Spatial distribution of the intensity field |p| 2 for the designed lens with/= 3X. (d) Transverse 
cross-section of intensity profile at x=63 cm for the designed lens. The black dash line refers to the intensity of the incident waves. 
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Figure 5 | Acoustic axicon for the non-diffracting Bessel beam, (a) Schematic diagram of the design of acoustic axicon. The green cone-like line with 
base angle = arctan(/z/DOF) is the desired equiphase surface for the acoustic axicon so that the phase shift at Pa should be proportional to the length of 
the red line P a Sa- Here, the half-height of the axicon h is selected to be 100 cm. (b) The theoretical continuous phase shift (red dots) and the discrete phase 
shift provided by the metasurface (blue squares) along the y direction, (c) Spatial distribution of the intensity field \pf for the designed axicon with 
DOF = 28/l. (d) Transverse cross-section of the intensity profile at x = 240 cm from the metasurfaces within the DOF. 



derivation of the generalized Snell's law. Considering the acoustic 
waves with normal incidence along the — x direction, due to the 
existing of discrete phase shifts along the y direction, the angle of 
the reflected wave 0 r (measured from the x direction) can be deduced 
by applying Fermat's principle 1 : 



, = arcsin 



where (j){y), dy represent the phase shift and the distance between two 
cross points along the y direction respectively, and k = 2n/X is the 
wave vector in air. Equation (1) implies that the reflected angle can be 
engineered freely by designing the suitable phase profile along the y 
direction. 

Suppose that the gradient of the phase profile <j>(y) on the metasur- 
face is constant along the y direction. From Eq. (1), it can be con- 
cluded that the reflected waves should remain a plane wave and be 
independent of the location y. To demonstrate this feature, one type 
of phase profile is introduced. The phase profile is chosen as 
d(j)/dy = n/16p, indicating that there are four identical units for each 
discrete phase shift of tc/4. The pressure field pattern p(x, y) is 



illustrated in Fig. 2(b). Different from normal surfaces, by which 
the reflected waves can only propagate along the opposite direction 
of the normal incidence, anomalous reflection can be clearly 
observed after imposing the discrete phase profile. The theoretical 
angles of reflection for this case should be 0 r = 37.77°, which can be 
deduced from Eq. (1). Excellent agreement can be observed by com- 
paring the theoretical value (shown as black arrows) and simulated 
pattern. Note that imperfect plane waves are obtained due to the 
mismatched pressure amplitude of the reflected acoustic waves along 
the metasurface as mentioned above. 

Moreover, we can design the specific distributions of <j> (y) = 
h sin 0 r -y to realize different wavefront manipulations. One interest- 
ing phenomenon is the conversion of a propagating wave into a 
surface wave. It can be verified from Eq. (1) that the reflected waves 
with 0 r = 90° could be obtained with the phase profile <j>(y) = ± hy. 
Different from the Ref. 4, where (j>(y) = h\y\ is utilized to realize the 
directional surface waves, the phase profile (j)(y) = — hy is employed 
to generate two surface beams propagating along they and —y direc- 
tions, respectively. For the sake of clarity, it is necessary to build a rule 
to convert a continuous phase profile to a discrete phase profile. 
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Throughout the paper, the rule is built as follows: for example, the 
units providing n/4 phase shift is utilized to displace the continuous 
phase range n/4 ± n/S at the specified position in they direction. The 
pressure field pattern p(x, y) of the reflected waves by the metasur- 
faces constructed from the rule is illustrated in Fig. 3. It is found that 
the surface waves are propagating near the metasurfaces along the ± 
y directions (black arrows) after the incident wave impinge on the 
metasurfaces. High efficiency conversion of propagating waves to 
surface waves is realized even though there is some diffraction upon 
the metasurfaces. Physically, the discrete phase shifts along the y 
direction provide extra momenta to compensate the momentum 
mismatch between the propagating wave and the surface wave on 
the metasurfaces, resulting in the high efficiency conversion. 

Next, we show that a flat lens and a flat axicon can be achieved 
successfully through acoustic metasurfaces. To design an acoustic flat 
lens, the hyperboloidal phase profile is employed on the metasur- 
faces, as shown in Fig. 4(a). For a given focal length/, the phase shift 
<j)(y) imposed at every point Pl(7) must satisfy the following 
equation 

*W = *-PE5>y (V7+P-/). (2) 

By applying the aforementioned rule, the desirable continuous phase 
profile (red dots) and the discrete phase profile (blue squares) pro- 
vided by the metasurfaces are plotted in Fig. 4(b). The physical struc- 
ture of the lens is readily constructed according to the phase profile. 
Spatial intensity distribution \p\ of the acoustic lens with / = 3X is 
shown in Fig. 4(c). To quantify the performance of the acoustic lens, 
the transverse cross -section intensity distribution at x = 63 cm is 
shown in Fig. 4(d). The intensity of pressure at the focal point is 
nearly 14.13 times larger than the one of the incident waves (shown 
as black dash line), which provides clear confirmation that excellent 
focusing effect can be realized by the proposed metasurfaces. 
Note that the hyperboloidal phase profile imposed along the y 
direction could produce a spherical wavefront even for nonpara- 
xial conditions, which leads to excellent focusing effect without 
aberrations 5 . 

To design a flat acoustic axicon with a base angle the conical 
phase profile should be employed so that the phase shift (j)(y) at every 
point Pa (y) along the y direction satisfy the following equation 



<Ky)=/cP A SA=yMsinjS (3) 

By applying the aforementioned rule, the theoretical continuous 
phase shift (red dots) and the discrete phase shifts (blue dots) 
provided by the acoustic metasurfaces are shown in Fig. 5(b). 
Spatial intensity distribution |p| 2 of the acoustic axicon and the 
transverse cross-section of the intensity profile at x = 240 cm 
from the metasurface are shown in Figs. 5(c) and 5(d), respect- 
ively. Due to the good performance of the proposed acoustic 
metasurfaces, it is not surprising to observe a non- diffracting 
Bessel beam propagating along the x direction with a relatively 
long distance. 

Discussion 

In summary, we have demonstrated that acoustic metasurfaces can 
be constructed by eight units of corrugated plates, which can tailor 
the reflected waves with discrete phase shifts covering the full In 
span. The generalized Snell's law can be extended to acoustic 
counterpart, which is otherwise difficult to realize with natural mate- 
rials. By employing linear phase profiles, anomalous reflections and 
the conversion of propagating waves into surface waves can be rea- 
lized. Furthermore, a planar acoustic lens and an acoustic axicon can 
be constructed by utilizing the hyperbolic and conical phase profiles. 
Our results should open the door for wavefront manipulation and 
engineering of acoustic waves. 



Methods 

Throughout the paper, Finite Element Method (FEM) based on commercial software 
COMSOL Multiphysics™ 4.3a is employed for the simulations. The materials applied 
in the simulations are air and stiff beams (sound hard boundaries). Plane wave 
radiation boundary condition is imposed on the incident boundaries and the periodic 
boundary condition are employed in the y direction to calculate the phase profile of 
the reflected wave for different length a (Fig. l(c-d)) and the pressure field patterns 
(Fig. 2(b)). Each strips of pressure pattern shown in Fig. 1(d) is calculated indepen- 
dently. For Figs. 3-5, perfectly matched layers (PMLs) and background pressure field 
is utilized to eliminate the reflected waves by the outer boundaries because of the 
plane wave radiation boundaries are not sufficient to absorb the waves with oblique 
incidence. 
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